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Abstract—It is known from both theory and numerical simula- 
tions that a current pulse suffers apparent attenuation as it propa- 
gates along a vertical perfect conductor of uniform, nonzero thick- 
ness (e.g., a cylinder) above perfectly conducting ground, excited 
at its bottom by a lumped source. The associated electromagnetic 
field structure is non-transverse electromagnetic (TEM), particu- 
larly near the source region. On the other hand, it has been shown 
analytically by Thottappillil et al. (2001, 2004) that no attenuation 
occurs and the electromagnetic field structure is pure transverse 
electromagnetic (TEM) if the conductor thickness and source size 
are assumed to be infinitesimal. The goal of this paper is to exam- 
ine the mechanism of current attenuation as it propagates along a 
nonzero thickness conductor, based on the scattering theory and 
on a nonuniform transmission line approximation. In applying the 
scattering theory, we decompose the “total” current in the con- 
ductor into two components that we refer to as the “incident” and 
“scattered” currents. The “incident” current serves as a reference 
(no attenuation), specified disregarding the interaction of resultant 
electric and magnetic fields with the conductor, while the “scat- 
tered” current, found here using the finite-difference time-domain 
(FDTD) method, can be viewed as a correction to account for that 
interaction. The scattered current modifies the incident current so 
that the resultant total current pulse appears attenuated. Thus, 
the current attenuation is likely to be due to field scattering that 
does not occur in the case of zero thickness conductor. The attenu- 
ation of the total current pulse is accompanied by the lengthening 
of its tail, such that the total charge transfer is independent of 
height. Approximation of the vertical conductor above ground by 
a nonuniform transmission line whose characteristic impedance 
increases with increasing height is shown to reasonably reproduce 
the current pulse attenuation predicted by the scattering theory. In 
this approximation, the apparent current attenuation with height 
can be attributed to waves reflected back to the source. The results 
have important implications for development and interpretation of 
lightning models. 


Index Terms—Current attenuation, lightning, TEM, transmis- 
sion line model, vertical conductor. 


I. INTRODUCTION 


HE METHOD of moments (MOM) (Harrington [1]) and 
T the finite-difference time-domain (FDTD) method (Yee 
[2]) for solving Maxwell’s equations have recently been ap- 
plied to finding the distribution of current along the light- 
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ning channel. Models based on this approach are referred to 
as electromagnetic models (e.g., Rakov and Uman [3]). In these 
models (e.g., Podgorski and Landt [4], Moini et al. [5], Kordi 
et al. [6], Baba and Rakov [7]), lightning return stroke chan- 
nel is represented by a vertical conductor of nonzero thickness. 
It has been shown that a current pulse attenuates as it prop- 
agates along such a conductor, even in the absence of ohmic 
losses, when it is excited at its bottom by a lumped voltage 
source (Bantin [8]; Kordi et al. [9], [10]; Baba and Ishii [11]) 
or by a lumped current source (Grcev ef al. [12], and Baba 
and Rakov [7]). For example, the current peak decreases from 
11 kA to 10 kA over the bottom 300 m of a vertical perfectly 
conducting wire of radius 5 cm in air above perfectly conducting 
ground (Kordi et al. [10] and Grcev et al. [12]), when a current 
expected for lightning subsequent return strokes (Nucci et al. 
[13]) is injected into the wire. The current attenuation predicted 
by numerical simulations is consistent with theory (e.g., Wu 
[14], Chen [15], Lee and Dudley [16]), as shown in Appendix II. 

On the other hand, it has been shown analytically (Thottap- 
pillil et al. [17], [18]) that no attenuation occurs and electro- 
magnetic field structure is pure TEM if the conductor thickness 
and source size are assumed to be infinitesimal. This case corre- 
sponds to the so-called transmission line (TL) model (Uman and 
McLain [19]), not to be confused with distributed-circuit (R-L- 
C) models (e.g., Rakov and Uman 1998 [3]). The TL model 
with return-stroke speed, v, equal to the speed of light c can 
be visualized, based on the transmission line theory, as a zero- 
radius vertical perfectly conducting wire and its image excited 
at their junction point by an infinitesimal source, as shown in 
Fig. 1(a). The image can be viewed as the return path for the ver- 
tical wire, and vice versa. An equivalent representation, based 
on the waveguide theory, is a zero-radius vertical perfectly con- 
ducting wire above perfectly conducting ground excited at its 
bottom by a spherical TEM wave (Kordi et al. [10]), as shown 
in Fig. 1(b). These two visualizations require a hypothetical 
zero-radius vertical wire, which cannot be realized in practice 
or in numerical techniques. The most general and rigorous rep- 
resentation of the TL model, based on the antenna theory, is 
a vertical phased current source array (Baba and Rakov [7]), 
shown in Fig. 1(c). Each current source is activated when an 
upward moving return stroke front reaches its altitude. This lat- 
ter representation is valid for any return stroke speed v and any 
channel radius. 

There appears to be no consensus regarding the reason (or 
reasons) for current attenuation predicted by the electromag- 
netic lightning models. According to Kordi et al. [9], [10], 
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Fig. 1. Visualization of the lightning return stroke TL model, [(z',t) = 


I(0,t — z'/v), with v = c. Semicircular solid lines and slanted arrows rep- 
resent an outward propagating spherical TEM wave front. (a) Visualization 
based on the classical transmission line theory: a vertical wire and its image 
excited by an infinitesimal source, valid only for v = c (v < c would imply 
a medium other than air) and a zero radius lightning channel (ZL = o0,C' = 0 
where L and C are the inductance and capacitance, each per unit length). A 
spherical TEM field structure is required for the classical transmission line the- 
ory to be applicable. (b) Visualization based on the waveguide theory: a vertical 
wire above ground excited at its bottom by a spherical TEM wave, valid only for 
v =c(v < c would imply a medium other than air) and a zero radius lightning 
channel. (c) Visualization based on the antenna theory: a vertical phased current 
source array above ground, valid for any lightning return stroke speed, v, and 
any lightning channel radius. This is the most general and rigorous representa- 
tion of the TL model. Note that for the case of v < c, the electromagnetic field 
structure is not TEM. 


the current attenuation and dispersion necessarily occur when 
a vertical nonzero radius wire is excited by a non-TEM- 
wave source. Baba and Rakov [7] simulated spherical TEM 
wave excitation using an appropriate phased current source ar- 
ray and observed current attenuation above the excitation re- 
gion. In view of the absence of ohmic losses, the observed 
current attenuation was attributed to “radiation losses” (e.g., 
Kordi et al. [9], [10]; Baba and Ishii [11]). Thottappillil and 
Uman [20] suggested that a vertical cylinder above ground can 
support an unattenuated current if this cylinder is excited at its 
base by a circular ensemble of TEM wave sources. On the other 
hand, Baba and Rakov [7] have pointed out that the charac- 
teristic impedance of a vertical cylinder above ground varies 
with height, which implies that pulses propagating along such a 
cylinder in general cannot preserve their amplitude and shape. 
Bermudez et al. [21] stated that the controversy regarding the 
validity of uniform transmission line representation for a ver- 
tical conductor above ground “is not settled yet.” One of the 
reasons for the contradictory views found in the literature is 
the lack of understanding of the mechanism of attenuation of 
current waves propagating along a perfectly conducting wire. 
In this paper, we will attempt to explain the mechanism of 
the current attenuation based on the scattering theory and on 
a nonuniform transmission line approximation. In applying the 
scattering theory, we decompose an attenuated “total” current 
pulse /,,, propagating on a vertical nonzero-thickness perfect 
conductor into an “incident” unattenuated current pulse, inc, 
and an induced or “scattered” current pulse, Jscat (evaluated 
as Itot — Linc). This is equivalent to finding a correction to the 
initially assumed reference (no attenuation) current, Ji,., in the 
conductor, this correction, Isca;, accounting for the interaction 
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Fig. 2. Conceptual picture to explain the mechanism of current attenuation 
along a vertical nonzero thickness perfect conductor above perfectly conducting 
ground. All currents are assumed to flow on the axis. An attenuated “total” 
current pulse J¢,¢ is separated into an “incident” unattenuated current pulse 
Tinc and a “scattered” current pulse Iscat- [inc generates an incident downward 
vertical electric field at a horizontal distance x from the axis (on the lateral 
surface of the cylinder). Iscat produces a scattered upward vertical electric field 
that cancels the incident downward vertical electric field on the surface of the 
cylinder, and modifies the incident current J;,.. The resultant current pulse, 
Ttot = Tine + Iscat, appears attenuated and its tail is lengthened as this pulse 
propagates along the conductor. 


of resultant electric and magnetic fields with the conductor. Each 
of the three currents is assumed to flow on the axis of the conduc- 
tor and is represented by an appropriate vertical phased current 
source array (see Fig. 2). Our goal is to examine the mechanism 
of current attenuation along a nonzero-thickness vertical per- 
fectly conducting wire; in other words, to show the reason why 
a nonzero-thickness vertical perfectly conducting wire (as op- 
posed to a zero-thickness wire) cannot support “unattenuated” 
current propagation. This is why we introduce “unattenuated” 
incident current jn. as the reference point. An incident unatten- 
uated current J;,- transporting positive charge upward gener- 
ates an incident downward vertical electric field. We will show 
that the corresponding scattered current pulse Ica, produces an 
upward-directed scattered electric field that, as expected, can- 
cels the incident downward vertical electric field on the surface 
of the perfect conductor. This result indicates that the current 
pulse attenuation occurs in order to satisfy the boundary condi- 
tion on the tangential component of electric field on the surface 
of a perfect conductor. Fig. 2 shows the conceptual picture of 
this mechanism. Note that placing the equivalent filamentary 
current on the axis and matching the field boundary conditions 
on the surface is equivalent to assuming the equivalent filamen- 
tary current on the surface and placing matching points on the 
axis (e.g., Balanis [22]). In the nonuniform transmission line 
approximation, the apparent current attenuation with height can 
be attributed to waves reflected back to the source. 

The structure of this paper is as follows. In Section II, us- 
ing Thottappillil et al.’s [17] analytical equation, we examine 
the incident vertical electric field, Ein., which can be viewed 
as generated by an incident unattenuated current pulse Jj,,. In 
Section III, using the FDTD method, we find the total current 
Ito, needed to satisfy the boundary condition on the electric 
field on the surface of a vertical conductor of nonzero thickness. 
We show that J;,4 is attenuated, and that its shape changes, so 
that no charge is deposited on the conductor. In Section IV, 
we find the scattered current J,.,; as the difference between 
the total and incident currents (Ito: — Jinc). We examine the 
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relation of [,.,4, which can be viewed as a correction to the ini- 
tially assumed unattenuated current, Jjn-, to the attenuation of 
the total current pulse and lengthening of its tail. In Section V, 
we examine the dependencies of current attenuation on conduc- 
tor thickness, frequency, and source length. In Section VI, we 
show that a nonuniform transmission line whose characteristic 
impedance increases with increasing height can reasonably re- 
produce the current pulse attenuation observed using the FDTD 
method in Section III. In Appendix I, we explain the reason for 
less attenuation in the case of a horizontal conductor. In Ap- 
pendix II, we compare total current waveforms calculated using 
the FDTD method with those calculated using an analytical 
equation derived by Chen [15] for an infinitely long, perfectly 
conducting cylinder excited by a delta-gap step voltage source. 
In Appendix III, we examine the electric field structure around a 
nonzero-thickness perfect conductor. In Appendix IV, we show 
the dependencies of magnitude of the longitudinal electric field, 
produced by a short current dipole, on horizontal distance from 
the current dipole and on pulse width, which are used in Sec- 
tion V to explain the dependencies of current attenuation on 
conductor thickness and on frequency. 


II. INCIDENT CURRENT (Jinc), INCIDENT E-FIELD (Finc): 
ANALYTICAL SOLUTION 


In the TL model, the longitudinal current I(z’, t) in a straight 
and vertical lightning channel at an arbitrary height z’ and time 
t is expressed as follows: 


I(z',t) = 1(0,t — 2’/v) (1) 


where v is the return-stroke speed and [(0, t) is the channel-base 
current. Typical values of v are one-third to two-thirds of c (e.g., 
Rakov [23]). In this paper, we will assume that v = c, as done 
by Thottappillil et al. [17], Kordi et al. [9], [10]; and Baba and 
Rakov [7], in order to simplify the analysis aimed at identifying 
the mechanism of current attenuation with height. 

Baba and Rakov [7] showed numerically that a vertical phased 
current source array activated so as to satisfy (1) with v = c pro- 
duced a spherical TEM wave [see Fig. 1(c)]. The corresponding 
analytical solution is found by Thottappillil et al. [17], with the 
electric field equation being given by 


1 
EB 6,t) = 
a(r,6,¢) 2Qregcer sind 


I(0,t —r/c), 
where r is the inclined distance from the channel origin to an 
observation point, and 6 is the angle between the channel and 
a straight line passing through both the channel origin and the 
observation point. Note that (2) accounts for the electrostatic, in- 
duction, and radiation field components (Thottappillil et a/. [17], 
Fig. 2). The r~! dependence of the total electric field is due to 
the assumption v = c (return stroke front moves at the same 
speed as electromagnetic signal). We will use a vertical, zero- 
thickness phased current source array to generate an incident 
spherical TEM wave. We will refer to the source current distri- 
bution [given by (1)] as the “incident” unattenuated current Iinc. 
This current is shown in Fig. 3. Although other Jin. are possi- 
ble for the same J;,,, our choice is the best for examining the 
mechanism of current attenuation with height. Then, in Section 


640 (2) 


523 


Current [ kA ] 


600 
Time [ns ] 


400 


Fig. 3. Incident current, J;,-, waveforms at different heights produced by a 
vertical phased current source array, each source generating a pulse having a 
magnitude of | kA and a half peak width of 67 ns. 
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Fig. 4. Phased current source array of zero thickness above perfectly conduct- 
ing ground, activated as specified by (1) with v = c. FE, (a, z,t) represents the 
incident vertical electric field Fine. 


TABLE I 
MAGNITUDES OF DOWNWARD VERTICAL ELECTRIC FIELD, Fine, INkV/m AT A 
RADIAL DISTANCE «x (@ = O MAY BE VIEWED AS x — 0) FROM THE AXIS OF A 
VERTICAL PHASED CURRENT SOURCE ARRAY AND AT A HEIGHT z ABOVE THE 
GROUND PLANE, CALCULATED USING (3). CURRENT MAGNITUDE IS SET AT 
1 kA. DOWNWARD FIELD IS DEFINED AS POSITIVE 


z,m 
x,m 0 0.1 1 10 100 
0 00 600 60 6.0 0.60 
0.1 600 420 60 6.0 0.60 
1.0 60 60 42 6.0 0.60 


III, we will consider a vertical nonzero-thickness conductor and 
find the distribution of current [,,.; along this conductor needed 
to satisfy the boundary condition on the electric field on its sur- 
face, using the FDTD method. The difference between /;.4 and 
Tine can be viewed as the scattered current [scat = Ltot — Linc. 

The vertical electric field component £, at a horizontal dis- 
tance x from the axis and at a height z above ground plane is 
given by 


E. (a, z,t) = Ey(r,0,t) cos (5 = 6) 


1 
= I(0,t — 
2m e9cr 03 ro) 


1 
= T(0,t-— Va224+22/c) 3 
QreqcV x2 + 22 ( fey) 
and illustrated in Fig. 4. 
Table I shows magnitudes of downward vertical electric field 
calculated using (3) along and around the vertical zero-thickness 
phased current source array, which represents an incident 


524 IEEE TRANSACTIONS ON ELECTROMAGNETIC COMPATIBILITY, VOL. 47, NO. 3, AUGUST 2005 


Cell size: 1x1x1m° 
> <40m 


Perfectly matched layers 
(10 m thick) 


Perfectly conducting vertical 
rectangular parallelepiped 
(2 mx 2 mx 300 m) 


300 m 


Air : €0, Ho 


Perfectly conducting plane 


~ 1-m high lumped current source 


Fig. 5. A vertical perfectly conducting parallelepiped having a cross sectional 
area of 2 x 2 m? in air ona perfectly conducting plane excited at the bottom by a 
1-m high current source. Current, [¢t, in this conductor corresponding to zero 
vertical electric field (E404 = 0) on the surface of conductor is found using 
the FDTD method. The working volume of 40 x 40 x 310 m? is surrounded by 
perfectly matched layers (PML) of thickness 10 m, except for the bottom, where 
the volume is limited by the perfectly conducting plane. 


unattenuated positive current pulse Ji,, having a magnitude of 
1 kA and propagating upward with v = c. Interestingly, the tan- 
gential (vertical) component of electric field at x = 0 is nonzero 
(approaches a nonzero value as x approaches zero, x — 0), im- 
plying that the phased current source array cannot be viewed as 
a vertical perfectly conducting cylinder of zero radius. In this 
case, the TEM field structure is maintained due to the fact that, 
although the vertical component of electric field at x = 0 m is 
nonzero, its radial component (the same as the theta component 
at x = 0 m) is infinitely large. The downward vertical electric 
field along a zero-thickness array (x = 0 m, z # 0) varies as 
inverse height (1/z). Away from the array (2 > 0 m), the verti- 
cal electric field decreases with increasing height more slowly 
than the inverse proportionality at relatively small heights (near 
the source region), and decreases approximately as the inverse 
height beyond that region. 

Note that for two parallel vertical phased current source ar- 
rays, the vertical electric field decays with height much faster 
than inverse height in the case of a single vertical phased current 
source array and vanishes within distances exceeding several 
times the spacing between the two parallel current source arrays 
(Appendix I-A). This observation has implications for the case 
of horizontal wires above ground, considered in Appendix I-B. 


III. TOTAL CURRENT (Jto¢), TOTAL E-FIELD (E tot): 
NUMERICAL SOLUTION 


In order to simulate a spherical TEM wave modified by the 
presence of nonzero-thickness conductor and find associated 
current, J,,4, in the conductor (see Fig. 2), we consider a verti- 
cal perfectly conducting parallelepiped having a cross sectional 
area of 2 x 2 m? in air on a perfectly conducting plane, excited 
at its bottom by a | m high current source (Fig. 5). The solution 
will be found using the FDTD method designed for rectangular 
geometry. Conclusions based on this solution are also applicable 
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Fig. 6. Specification of 1 m high current source at the bottom of the 2 x 2 x 
300 m? vertical parallelepiped shown in Fig. 5. 
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Fig. 7. Total current waveforms, Jto¢, at different heights calculated using the 
FDTD method for a vertical perfectly conducting parallelepiped having a cross 
sectional area of 2 x 2 m? in air on a perfectly conducting plane excited at its 
bottom by a 1 m high lumped current source. This current source produces a 
Gaussian pulse having a magnitude of 1 kA and a half peak width of 67 ns. 


to the cylindrical geometry discussed above. The total current 
Jot minus the incident current J;,,. will yield the scattered cur- 
rent Iscat. AS noted above, this procedure can be viewed as 
finding a correction, Jsca¢, to the initially assumed unattenuated 
current, J;,-, needed to account for the reaction of the conductor. 
Alternatively, sca, could be obtained setting Fgcat = —Einc on 
the surface of the conductor and solving for J, on its axis. 

The working volume of 40 x 40 x 310 m3, shown in Fig. 5, is 
divided into 1 x 1 x 1 m° cubic cells. Perfectly matched layers 
(PML) (Berenger 1994 [24]) (absorbing boundaries) are set at 
the top and sides of the volume in order to avoid reflections there. 
At the bottom, the volume is limited by a perfectly conducting 
plane that is simulated by forcing the tangential components of 
electric field at the boundary to be zero. This latter method is also 
used to simulate vertical perfectly conducting parallelepiped. 

The current source is represented, as illustrated in 
Fig. 6, by specifying 12 1-m-long magnetic-field elements, 
H,,1(t), Hyi(t), and so on (that is, 3 x 3 m? magnetic-field 
loop), surrounding the bottom 2x 2x1 m?® section of the 
vertical parallelepiped. 

Fig. 7 shows current waveforms at different heights calcu- 
lated using the FDTD method with a time increment of 1.25 
ns for the vertical perfectly conducting parallelepiped shown in 
Fig. 5. The current source produces a Gaussian pulse having a 
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Fig. 8. Time integrals of current at heights 0, 10, and 50 m calculated using the 


FDTD method for a 1200 m long vertical perfectly conducting parallelepiped 
having a cross sectional area of 2 x 2 m? in air on a perfectly conducting plane 
excited at its bottom by a 1 m high lumped current source, which produces a 
Gaussian pulse having a magnitude of 1 kA and a half peak width of 67 ns. 


magnitude of | kA and a half-peak width of 67 ns. The current 
pulse propagates upward with attenuation, which is most signif- 
icant within the bottom 20 m, while the risetime of the current 
pulse is almost constant. Within the bottom 20 m, the average 
propagation speed is 0.9c, evaluated tracking the pulse peak, or 
0.95c, evaluated tracking the pulse rising edge that is defined as 
1% level of the corresponding peak on the rising edge, and is 
essentially equal to c at larger heights. 

As seen in Fig. 7, the current pulse attenuation is accompa- 
nied by the lengthening of its tail. Fig. 8 shows the time integral 
of current at heights 0, 10, and 50 m calculated using the FDTD 
method for a 1200-m long vertical perfectly conducting par- 
allelepiped having a cross-sectional area of 2 x 2 m? excited 
at its bottom by a 1-m high lumped current source. It is clear 
from Fig. 8 that, as time approaches infinity, the charge trans- 
ferred through each section of the conductor tends to be equal 
to the charge supplied by the source. In other words, as time 
approaches infinity, charge on the conductor approaches zero. 
This result possibly suggests that no energy is lost from the sys- 
tem (composed of an infinitely long vertical conductor above an 
infinitely large ground plane). If so, then, the apparent current 
attenuation, often viewed as being due to “radiation losses,” can 
be attributed to the redistribution of energy between electric and 
magnetic fields that is needed to satisfy the boundary conditions 
on the surface of perfect conductor. An apparently equivalent in- 
terpretation in terms of the nonuniform transmission-line theory 
is presented in Section VI. 


IV. SCATTERED CURRENT (Igcat,), SCATTERED E-FIELD 
(Escat) : I scat = Leot _ Linc; Escat = — Fine 


In this Section, we find the scattered current as Igca, = 
Itot — Tinc, which is responsible for the apparent attenuation 
of the total current /,,¢ with height. The total current wave- 
forms J,., at different heights, calculated using the FDTD 
method for a vertical perfectly conducting parallelepiped hav- 
ing a cross sectional area of 2 x 2 m? in air, are shown in 
Fig. 7. The corresponding incident current waveforms Inc, 
given by (1) for v = c, are shown in Fig. 3. Fig. 9 shows wave- 
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Fig.9. Scattered current waveforms, scat, at different heights obtained as the 
difference between corresponding waveforms shown in Figs. 7 and 3 (Iscat = 
iT tot — I; inc ): 
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Fig. 10. Waveforms of vertical electric field at a horizontal distance x = 1m 
from the conductor axis and at heights z = 10, 20, and 50 m calculated using 
appropriate phased current source arrays. Downward electric field is defined 
as positive. The positive unipolar waveforms are calculated using J;,. shown 
in Fig. 3 (these waveforms can also be obtained using (3)), and the negative 
unipolar waveforms are calculated using [scat shown in Fig. 9. Each source of 
the phased current source array is represented by a 1 x 1 m? magnetic field loop 
surrounding the conductor axis. 


forms of Iscat = Ltot — Tine (the difference between the wave- 
forms shown in Figs. 7 and 3). The scattered current [5.44 has 
a bipolar waveshape with the initial half cycle having a polarity 
that is opposite to that of the incident current J;,.. Thus, [scat is 
responsible for the apparent attenuation of /,._ with height. The 
second half cycle of Ica has the same polarity as Jj,- and causes 
the lengthening of the total current pulse tail. Since [4 deposits 
no charge on the vertical perfect conductor (see Fig. 8), the time 
integral of J,.a4 1S equal to zero (as time approaches infinity). 
Fig. 10 shows waveforms of the vertical electric field at a hor- 
izontal distance x = 1 m from the conductor axis and at heights 
z = 10, 20, and 50 m in air due to Jin (Fig. 3) and [scat (Fig. 9). 
The scattered electric fields E'scat (negative unipolar waveforms 
in Fig. 10) were calculated using the scattered current distri- 
bution Jscaz (Fig. 9) as the outputs of a vertical phased current 
source arrays (see Fig. 2). The incident electric fields Ein. (pos- 
itive unipolar waveforms in Fig. 10) were calculated similarly 
using the incident current distribution Ji,, (Fig. 3) as the outputs 
of a vertical phased current source arrays (or, equivalently, using 
(3)). Each current source was represented by specifying four 
1-m long magnetic-field elements (1 x 1 m? magnetic-field 
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Fig. 11. Current waveforms at different heights calculated using the FDTD 


method for a vertical perfectly conducting parallelepiped having a cross sec- 
tional area of 2 x 2 m? in air on a perfectly conducting plane excited at its 
bottom by a | m high lumped current source. This current source produces a 
current pulse representative of lightning subsequent return strokes. See Nucci 
et al. [13]. 


loop) surrounding the conductor axis, as described above. The 
scattered current J,.,4, which is responsible for the reduction 
of Io, pulse magnitude, produces a scattered upward vertical 
electric field that, as expected, cancels the incident downward 
vertical electric field 1 m away from the conductor axis. Thus 
the current attenuation is likely to be due to field scattering that 
does not occur in the case of zero-thickness conductor. As noted 
above, Iscaz is also responsible for the lengthened wavetail of 
the total current Jyot = Tine + Iscat. It follows from Table I and 
(3) that a current pulse attenuates significantly near the source 
region, and that attenuation reduces with increasing height. 

Both incident and scattered fields contribute to the overall 
electromagnetic field structure. Near the vertical conductor, the 
vertical electric field due to [scat decays with increasing hori- 
zontal distance x considerably faster than that due to Jinc. As 
a result, the TEM field structure due to Jin, is disturbed in the 
vicinity of the conductor. At larger distances, the fields due to 
Tine and Isca, decay at about the same rate, so that the overall 
electromagnetic field structure becomes more or less TEM (see 
Appendix II). 

Note that the same solution for [s¢,4 would be obtained if we 
placed a large number of phased current source arrays, whose 
total current equals to J;,-, on the lateral surface of the cylinder 
and required a zero total vertical electric field on the axis of 
the cylinder. This formulation is apparently equivalent to the 
vertical cylinder excited at its base by a circular ensemble of 
TEM wave sources considered by Thottappillil and Uman [20]. 


V. DEPENDENCES OF CURRENT ATTENUATION ON THE SOURCE 
LENGTH, CONDUCTOR THICKNESS, AND FREQUENCY 


Besides the base case, 1-m high current source, 2 x 2 m? 
cross-sectional area parallelepiped, 67 ns half-peak width Gaus- 
sian pulse, presented in Section III, we performed calculations 
also for the following conditions: 1) a 2-m high current source; 
2) a vertical perfectly conducting parallelepiped having a 1 x 
1 m? cross sectional area; and 3) a 33 ns half-peak width 
Gaussian pulse. Table IT shows the dependency of current atten- 
uation on the source height, Table III on the thickness of vertical 


TABLE II 
MAGNITUDES OF CURRENT PULSES, J++, AT DIFFERENT HEIGHTS (IN PERCENT 
OF THE MAGNITUDE AT 2/ = 0) CALCULATED USING THE FDTD METHOD FOR 
A VERTICAL PERFECTLY CONDUCTING PARALLELEPIPED HAVING A CROSS 
SECTIONAL AREA OF 2 x 2 m2. THE VERTICAL CONDUCTOR IS EXCITED AT ITS 
BOTTOM BY A | OR 2 m HIGH LUMPED CURRENT SOURCE, WHICH PRODUCES 
A GAUSSIAN PULSE HAVING A HALF-PEAK WIDTH OF 67 ns 


Source height, z’,m 

m 0 10 20 50 100 
1 100 74 68 61 56 
2 100 79 72 64 59 


TABLE III 
MAGNITUDES OF CURRENT PULSES, J+o4, AT DIFFERENT HEIGHTS (IN PERCENT 
OF THE MAGNITUDE AT 2/ = 0) CALCULATED USING THE FDTD METHOD FOR 
VERTICAL PERFECTLY CONDUCTING PARALLELEPIPEDS HAVING CROSS 
SECTIONAL AREAS OF 2 x 2 m? AND 1 x 1 m? ON A PERFECTLY CONDUCTING 
PLANE. THE VERTICAL CONDUCTOR IS EXCITED AT ITS BOTTOM BY Alm 
HIGH LUMPED CURRENT SOURCE, WHICH PRODUCES A GAUSSIAN PULSE 
HAVING A HALF-PEAK WIDTH OF 67 ns 


Cross-sectional 


area 0 10 20 50 


100 
100 79 73 66 61 
100 74 68 61 56 


TABLE IV 
MAGNITUDES OF CURRENT PULSES, J¢o+, AT DIFFERENT HEIGHTS (IN PERCENT 
OF THE MAGNITUDE AT 2/ = 0) CALCULATED USING THE FDTD METHOD FOR 
A VERTICAL PERFECTLY CONDUCTING PARALLELEPIPED HAVING A CROSS 
SECTIONAL AREAS OF 2 x 2 m?. THE VERTICAL CONDUCTOR IS EXCITED AT 
ITs BOTTOM BY A | m HIGH LUMPED CURRENT SOURCE, WHICH PRODUCES A 
GAUSSIAN PULSE HAVING A HALF PEAK WIDTH OF 33 OR 67 ns 


Half-peak z’,m 
width, ns 0 10 20 50 100 
33 ns 100 66 59 52 47 


67 ns 100 74 68 61 56 


conductor, and Table IV on the pulse width or the frequency 
content of excitation. It is clear from Tables II to IV that the cur- 
rent attenuation becomes more pronounced as the source height 
decreases, the thickness of vertical conductor increases, and the 
pulse width decreases (higher frequency content increases). 

As seen in Fig. 7, current attenuation is most pronounced in 
the bottom part of vertical conductor (near ground). Since there 
is an imposed current distribution (no field boundary condition 
to satisfy) over the source height, there is less pronounced 
attenuation in the case of a longer source at the bottom of the 
conductor. 

As seen in Table I, the magnitude of incident downward 
vertical electric field generated by Ji, depends weakly on the 
horizontal distance from the conductor axis when the horizontal 
distance is much smaller than the height above ground. On the 
other hand, the magnitude of an upward vertical electric field 
produced by a single-cycle bipolar current pulse, which roughly 
represents J,.,, and has a negative initial half cycle, along 
a short dipole, decreases rapidly as the horizontal distance 
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increases (see Appendix IV). Therefore, a larger J. 1s needed 
to compensate the incident downward vertical electric field 
on the conductor surface as the conductor radius becomes 
larger. This is the reason why current pulse attenuation is more 
pronounced for a larger thickness conductor. 

It is clear from (3) that the magnitude of the downward ver- 
tical electric field produced by J;,. along a vertical conductor 
is independent of frequency if v = c (Thottappillil et al. [17]). 
On the other hand, the magnitude of the vertical electric field 
produced by a short current dipole becomes smaller as the fre- 
quency gets higher (see Appendix IV). Within a few meters from 
the short dipole, the static component of electric field, which is 
related to the time integral of current, is dominant. Therefore, the 
electric field produced by the dipole at such horizontal distances 
(representing the thickness of the conductor) becomes smaller 
as the pulse width gets smaller (its higher frequency content 
increases) if the magnitude remains the same. As a result, in 
order to compensate the incident electric field on the conductor 
surface, the magnitude of J;ca, needs to be larger as the pulse 
width gets smaller (higher frequency content gets larger). This 
is the reason why a more narrow current pulse exhibits a more 
pronounced attenuation. 

The risetime of the Gaussian current pulse in Fig. 7 becomes 
only about 10% longer within the bottom 20 m or so than that of 
the injected waveform and remains essentially the same beyond 
that region. On the other hand, if we inject a current pulse 
representative of lightning subsequent return strokes (Nucci 
et al. [13]), which is characterized by a broader frequency spec- 
trum than the 67-ns wide Gaussian pulse, into the vertical per- 
fectly conducting parallelepiped shown in Fig. 5, the pulse rise- 
time becomes appreciably longer as the current pulse propagates 
upward (see Fig. 11). This is because higher frequency compo- 
nents are subject to considerably stronger attenuation than lower 
frequency components. 


VI. NONUNIFORM TRANSMISSION-LINE APPROXIMATION 


Now we will consider a nonuniform transmission line ap- 
proximation (e.g., Wagner and Hileman [25]; Menemenlis and 
Chun [26]) and show that it is consistent with our scattering 
theory analysis given in Sections I-IV. In this approximation, 
the apparent attenuation of current with height can be attributed 
to waves reflected back to the source. The effect is more pro- 
nounced near the bottom of the vertical conductor. 

Fig. 12 shows a lossless nonuniform transmission line ex- 
cited by a current source of zero length. Although Haase and 
Nitsch [27], [28] and Haase et al. [29] have recently derived 
telegrapher’s equations for nonuniform transmission lines in the 
frequency domain, we will use a simplified approach assuming 
that the vertical perfectly conducting wire above ground can 
be represented by a lossless transmission line whose character- 
istic impedance increases with height as 60 cosh _'(z'/r) (e.g., 
Cheng [30]). This expression gives the characteristic impedance 
in ohm of a horizontal perfectly conducting wire of radius r 
(not to be confused with the radial distance in (2) and (3) and 
in Fig. 4) at a height z’ above perfectly conducting ground. 
However, such an approach is often used in lightning mod- 
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Fig. 12. A lossless nonuniform transmission line, representing a vertical con- 
ductor above ground or, equivalently, a vertical conductor and its image, excited 
by a current source. The characteristic impedance Z, is assumed to vary with 
height z’ as 60 cosh~+ (z'/r). The radius of the vertical conductor r is set to 1 
m. The current source produces a Gaussian pulse having a magnitude of 1 kA 
and a half-peak width of 67 ns. 
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Fig. 13. Current waveforms at different heights calculated for the nonuniform 
transmission line shown in Fig. 12, to be compared with Fig. 7. Values of Z,. at 
10, 30, 50, and 100 m are 180, 246, 276, and 318 Q, respectively. The half peak 
width of a Gaussian pulse injected at the botom is 67 ns. 


eling (e.g., Bazelyan et al. 1978 [31]), since there is no exact 
equation for the characteristic impedance for a vertical nonzero- 
radius conductor above ground. In this analysis, the nonuniform 
transmission line shown in Fig. 12 is divided into 3 m uni- 
form transmission line segments whose constant characteristic 
impedances are given by 60cosh '(z'/r). The current source 
was assumed to produce a Gaussian pulse having a magnitude 
of | kA and a half-peak width of 67 ns (the same as in Sections 
II-IV). The travel time of current waves along a 3-m long seg- 
ment is 10 ns. A traveling wave analysis method, the so-called 
Bergeron method (Dommel [32]), solving standard telegrapher’s 
equations, was used to compute currents along such a dis- 
cretized nonuniform transmission line. The results are shown in 
Fig.13. The treatment of reflections in the Bergeron method 
is qualitatively illustrated in Fig. 14. Note that the cosh” 
expression for the characteristic impedance reduces to 
60In(2z'/r) when 2! >>r (cosh7!# = In[a + (a? — 1)/?] 
for x > 1), which yields almost the same current profile as 
that shown in Fig. 13. 

From comparison of Figs. 7 and 13, one can see that the 
lossless nonuniform transmission line shown in Fig. 12 can 
reasonably reproduce the current pulse attenuation and the 
lengthened wavetail, obtained from the full-wave analysis us- 
ing the FDTD method for the 2x2 m? vertical perfectly 
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Fig. 14. A transmission line composed of two lossless uniform transmission 
line sections connected in series, excited by a current source. The characteristic 
impedance of the section above junction point A is higher than that of the 
section below A. When an incident upward propagating positive current wave 
(#1) encounters junction A, a downward negative (reflected) current wave (#2) 
is produced there. When the downward negative reflected current wave (#2) 
encounters the current source (infinitely large impedance), an upward positive 
(reflected) current wave (#3) is produced. 


conducting parallelepiped above ground (see Section II). The 
observed small differences in current decay seen in Figs. 7 and 
13 are likely to be due to different cross sections used in the 
FDTD (2 x 2 m?) and nonuniform transmission line (1 m ra- 
dius circle) simulations. In the case of a shorter injected current 
pulse (half-peak width of 33 ns), the waveforms of a current 
pulse at different heights calculated using the nonuniform trans- 
mission line shown in Fig. 12 are in good agreement with those 
calculated using the FDTD method for the vertical perfectly 
conducting parallelepiped (neither of these waveforms is shown 
in this paper). 

We now qualitatively discuss the current pulse attenuation 
and lengthening of its tail, predicted by the nonuniform trans- 
mission line approach. Fig. 14 shows the bottom two sections of 
discretized nonuniform transmission line presented in Fig. 12. 
When an incident upward propagating positive current wave 
(#1) encounters junction A, a downward negative (reflected) 
current wave (#2) is produced there. The reflected current wave 
is negative because the characteristic impedance above A is 
greater than below A. The apparent current attenuation with 
height along the nonuniform transmission line (Fig. 13) can be 
attributed to such waves reflected back to the current source. 
When the downward negative reflected current wave (#2) en- 
counters the current source (infinitely large impedance), an up- 
ward positive (reflected) current wave (#3) is produced. This 
latter reflection is apparently responsible for the lengthened tail 
of the total current waveform (Fig. 13). 

The time variations of the time integral of current pulse at 
each height for this nonuniform transmission line, which are 
not shown in this paper, are essentially identical to those cal- 
culated using the FDTD method for the vertical perfectly con- 
ducting parallelepiped (Fig. 8). Note that an attenuated current 
propagating along a vertical conductor above ground produces a 
non-TEM wave (Baba and Rakov [7]). Therefore, a similarly at- 
tenuated current pulse, which propagates along the nonuniform 
transmission line, also produces a non-TEM wave. 

Since the Bergeron method employs standard telegrapher’s 
equations for each segment of the wire, it apparently does not 
account for electromagnetic radiation. This fact and the simi- 


larity of the results obtained using the nonuniform transmission 
line approximation and scattering theory appear to support our 
speculation (see Section III) that the current attenuation can be 
explained without invoking the not-well-defined concept of "ra- 
diation losses." As noted above, in the nonuniform transmission 
line approximation, the apparent current attenuation with height 
can be attributed to waves reflected back to the source. 


VII. SUMMARY 
A. Attenuation Mechanism in Terms of the Scattering Theory 


The mechanism of apparent attenuation of current pulses as 
they propagate along a vertical nonzero-thickness perfect con- 
ductor above perfectly conducting ground can be visualized as 
follows. A reference (no interaction with the conductor, no at- 
tenuation) positive current pulse Jj, propagating upward gen- 
erates an incident spherical TEM wave, with vertical electric 
field on the surface of the conductor being directed downward. 
Cancellation of this field, as required by the boundary condition 
on the tangential electric field on the surface of a perfect con- 
ductor, gives rise to a “scattered” current Jsca¢. This scattered 
current Ica, modifies Ji,., so that the resultant total current 
pulse J/,,.¢ appears attenuated as it propagates along the vertical 
conductor. Thus it appears that the field scattering (not present 
in the case of zero-thickness conductor) causes current attenu- 
ation with height. The attenuation of the total current pulse is 
accompanied by the lengthening of its tail, such that the total 
charge transfer is independent of height. The electromagnetic 
field structure associated with an attenuated current distribution 
along a vertical conductor above ground is non-TEM. 


B. Attenuation Mechanism in Terms of the Nonuniform 
Transmission Line Theory 


A nonuniform transmission line whose characteristic 
impedance increases with increasing height is shown to reason- 
ably reproduce the current pulse attenuation and lengthening of 
its tail predicted by the scattering theory. In this representation, 
the apparent current attenuation with height can be attributed to 
waves reflected back to the source. 


C. Dependence on Height 


The magnitude of the incident downward vertical electric 
field on the conductor surface produced by Iinc is approximately 
inversely proportional to height if the height is much larger than 
the conductor thickness. This is the reason for less pronounced 
current attenuation with increasing height. 


D. Dependence on Pulse Width 


The magnitude of the incident vertical electric field on the 
surface of a perfect conductor, produced by Jinc, is independent 
of frequency if v = c. On the other hand, the magnitude of the 
scattered vertical electric field is related to the time integral of 
Iscat- Therefore, if the pulse width of incident vertical electric 
field on the conductor surface, to be cancelled by the scattered 
vertical electric field, gets narrower while its magnitude is the 
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I(0,t-z’/c) 


- (0,t-z/c) 


Perfectly conducting ground 


Fig. 15. Two parallel vertical phased current source arrays of zero thickness 
above perfectly conducting ground. The arrays simulate incident unattenuated 
opposite polarity current pulses having an amplitude of 1 kA propagating upward 
with v = c. The spacing between the arrays is d = 1 m. 


same, a larger Iscaz is needed. This is the reason why more 
narrow current pulses exhibit a stronger attenuation. 


E. Implications for Lightning Models 


Our finding that a current pulse propagating along a verti- 
cal nonzero-thickness conductor attenuates significantly near 
the ground has important implications for lightning modeling, 
since this is the region of primary interest when the generation of 
lightning fields at early times is considered. Lightning return- 
stroke models with an imposed current distribution along the 
entire channel represent a useful engineering tool, but should be 
used with caution in studying lightning physical processes. A 
vertical conductor above ground, representing a lightning chan- 
nel or a tall strike object, generally cannot support unattenuated 
propagation of current waves, even if all losses are neglected. 


APPENDIX I 


A. Incident E-Field for Two Parallel Vertical Phased Current 
Source Arrays: Analytical Solution 


Fig. 15 shows two parallel vertical phased current source 
arrays of zero thickness on perfectly conducting ground. The 
arrays represent incident unattenuated opposite polarity current 
pulses having the same amplitude and propagating upward with 
v=C. 

Table V shows magnitudes of downward vertical electric field 
along the axis of the left array, calculated using (3). Note that 
retardation of wave propagation from the right array is neglected 
in this calculation. In Table V, magnitudes of downward vertical 
electric field in the case of a single array are also shown for 
reference. It follows from Table V that the vertical electric field 
(to be compensated by the scattered field) decreases much faster 
in the case of two arrays than in the case of a single array and 
vanishes within heights exceeding several times the spacing 
between the two arrays. This result indicates that a current pulse 
on a horizontal conductor above perfectly conducting ground 
would attenuate less than that on a single vertical conductor 
above ground. 
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TABLE V 
MAGNITUDES OF DOWNWARD VERTICAL ELECTRIC FIELD IN kV/m AT A 
HEIGHT z ALONG THE AXIS (2 = 0 OR 2 — 0) OF THE LEFT PHASED CURRENT 
SOURCE ARRAY ABOVE GROUND PLANE SHOWN IN Fic. 15, CALCULATED 
USING (3). RETARDATION OF WAVE PROPAGATION FROM THE RIGHT PHASED 
CURRENT SOURCE ARRAY IS NEGLECTED. MAGNITUDES OF DOWNWARD 
VERTICAL ELECTRIC FIELD IN THE CASE OF A SINGLE PHASED SOURCE ARRAY 
ARE ADDED FOR REFERENCE. DOWNWARD FIELD IS DEFINED AS POSITIVE 


Configuration z,m 

0.1 1 10 ~=100 
Two parallel vertical arrays | 540 18 0.03 0.00 
Single vertical array 600 60 6.0 0.60 


E-fields due to incident currents 
a 300 m 


. m x 1 m horizontal conductor 


1-m long horizontal current source 


300 m 
< 


Perfectly conducting ground 


(a) 
300 m__ E-fields due to incident currents 


> 300 m 


1 mx 1m horizontal conductor 


5.5m 
‘canoes 5-element vertical in-phase 
current source array 


Perfectly conducting ground 
(b) 


Fig. 16. A horizontal conductor having a cross sectional area of 1 x 1 m? 
placed 5.5 m above perfectly conducting ground, to be analyzed using the 
FDTD method: (a) A 1 m long current source is inserted in the middle of 600 
m long horizontal wire. (b) A 5 m vertical in-phase current source array is 
connected to the middle of a 600 m long horizontal wire. Each current source 
produces a Gaussian pulse having a magnitude of | kA and a half-peak width 
of 67 ns. The white arrows indicate the directions of longitudinal electric fields 
generated by incident unattenuated current pulses. 


B. Total Current for Horizontal Configurations: 
Numerical Solution 


Fig. 16 shows a horizontal perfectly conducting paral- 
lelepiped having a cross sectional area of 1 x 1 m? placed 
5.5 m above perfectly conducting ground, to be analyzed us- 
ing the FDTD method. The working volume for this analysis 
is 40 x 620 x 40 m? and is surrounded by perfectly matched 
layers (PML) of thickness 10 m, except for the bottom, 
where the volume is limited by the perfectly conducting 
plane. The overall rectangular volume shown in Fig. 16 
is divided into 1x 1x1 m® cubic cells. The white ar- 
rows in Fig. 16 indicate the directions of the longitudinal 
electric fields generated by the incident unattenuated current 
pulses on horizontal wires. 

In Fig. 16(a), which is quite similar to the configuration of 
Fig. 15, both right-directed positive and left-directed negative 
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Fig. 17. Current waveforms at different distances from the current injection 
point calculated using the FDTD method for a horizontal conductor placed 
5.5 m above perfectly conducting ground: (a) For the configuration of Fig. 16(a). 
(b) For that of Fig. 16(b). 


(incident unattenuated) current pulses will produce left-directed 
incident electric fields. In Fig. 16(b), the longitudinal incident 
electric field near the current injection point would be much 
less than that in Fig. 16(a). Indeed, the incident right- and left- 
directed positive current pulses will produce left- and right- 
directed incident electric fields, respectively, which will tend to 
cancel each other. 

Fig. 17 shows current waveforms at different distances from 
the current injection point calculated using the FDTD method 
for each of the two cases presented in Fig. 16. The FDTD 
simulation results are consistent with what we expected from 
the examination of longitudinal incident electric fields generated 
by the incident current pulses. This indicates that longitudinal 
electric fields generated by the incident current pulses are closely 
related to current attenuation. Note that, even in Fig. 17(a), a 
current pulse decays only within the first 20 m or so, and it 
propagates with little or no attenuation beyond this region, where 
the incident longitudinal electric field essentially vanishes. This 
is generally referred to as a “quasi-TEM” field structure. 


APPENDIX II 


COMPARISON OF FDTD SIMULATION WITH AN ANALYTICAL 
SOLUTION 


Chen [15] has derived an approximate analytical equation for 
the transient current J (z’, t) along an infinitely long, perfectly 
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Fig. 18. Current waveforms at different heights, calculated using the FDTD 
method for a vertical perfectly conducting parallelepiped having a cross sec- 
tional area of 2 x 2 m? above a perfectly conducting plane excited at its bottom 
by a | m high lumped voltage source and those calculated using (II1) for a ver- 
tical perfectly conducting cylinder of 1 m radius above a perfectly conducting 
plane excited at its bottom by a delta-gap voltage source. In both cases, the 
voltage sources produce a ramp front wave having a magnitude of | MV anda 
risetime of 100 ns. 


conducting cylinder in air, excited in the middle by a delta gap 
step voltage, V, source. This equation is reproduced below. 


He = see ( eee ) 
o 2 In(Vc24? — 2? /a) 
where 19 is the free space impedance (12077), In is the natural 
logarithm, and a is the radius of the cylinder. Chen’s [15] 
approximate expression yields quite accurate results that are 
almost identical to those given by the exact formula (Wu [14]). 
If we apply (1) to a vertical cylinder placed on a perfectly 
conducting plane and excited at its bottom, we have only to 
multiply the magnitude of resultant current by 2, in order to 
account for the image source. 

Fig. 18 shows current waveforms at different heights, 2’, cal- 
culated using the FDTD method for a vertical perfectly con- 
ducting parallelepiped having a cross-sectional area of 2 x 2 m? 
on a perfectly conducting plane excited at its bottom by a | m 
high lumped voltage source and those calculated using (II1) for 
a vertical perfectly conducting cylinder of | m radius excited at 
its bottom by a delta-gap voltage source. In both cases the volt- 
age sources are assumed to produce a ramp-front wave having 
a magnitude of 1 MV and arisetime of 100 ns. Since (II1) is the 
solution for a step voltage excitation, we obtained the response 
to our ramp-front voltage wave using numerical convolution. 
The waveforms calculated using the FDTD method agree rea- 
sonably well with those calculated using (II1). Note that the 
small difference in current waveforms at z’ = 0 m is probably 
due to the difference in source size: 1 min the FDTD simulation 
versus zero in (II1). It is also worth noting that as the radius, 
a, of the cylinder approaches zero, current predicted by (II1) 
approaches zero as well, unless voltage, V, is infinitely large. 


(II1) 


APPENDIX III 


E-FIELD STRUCTURE AROUND A VERTICAL 
NONZERO-THICKNESS PERFECT CONDUCTOR 


We now examine the electromagnetic field structure in air sur- 
rounding the 2 x 2 x 300 m? or 1 x 1 x 300 m® vertical perfect 
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TABLE VI 
MAGNITUDES OF ELECTRIC FIELD (IN PERCENT OF THE MAGNITUDE OF THE 
CORRESPONDING ELECTRIC FIELD OF SPHERICAL TEM WAVE) AT DIFFERENT 
POINTS (7, 9) (SEE FIG. 4) AROUND THE 2 x 2 x 300 m? VERTICAL 
CONDUCTOR (FIG. 5), CALCULATED USING THE FDTD METHOD. THE 
VERTICAL CONDUCTOR Is EXCITED AT ITS BOTTOM BY A | m HIGH LUMPED 
CURRENT SOURCE, WHICH PRODUCES A GAUSSIAN PULSE HAVING A HALF 
PEAK WIDTH OF 67 ns 


O= n/l4 O= n/2 
rm vertical horizontal vertical 
14 66 95 83 
35 63 76 2B 
71 63 70 71 
141 62 66 70 
TABLE VII 
MAGNITUDES OF ELECTRIC FIELD (IN PERCENT OF THE MAGNITUDE OF THE 


CORRESPONDING ELECTRIC FIELD OF A SPHERICAL TEM WAVE) AT 
DIFFERENT POINTS (7, 9) (SEE FIG. 4) AROUND THE 1 x 1 x 300 m®? VERTICAL 
CONDUCTOR (FIG. 5), CALCULATED USING THE FDTD METHOD. THE 
VERTICAL CONDUCTOR IS EXCITED AT ITS BOTTOM BY A | m HIGH LUMPED 
CURRENT SOURCE, WHICH PRODUCES A GAUSSIAN PULSE HAVING A HALF 
PEAK WIDTH OF 67 ns 


O= m/4 O= an/2 

rm vertical horizontal vertical 
14 74 100 90 
35 70 82 80 
71 69 76 77 
141 69 72 76 


conductor (Fig. 5). Tables VI and VII show the magnitudes of 
vertical and horizontal electric fields at different points around 
these vertical conductors, calculated using the FDTD method. 
Note that, in calculations presented in this Appendix, the spac- 
ing between the vertical conductor axis and the right-hand side 
absorbing boundary wall in Fig. 5 was moved from 20 to 150 m. 
In a spherical TEM wave, vertical and horizontal electric fields 
at 0 = 7/4, and vertical electric field at 9 = 7/2 should be the 
same. One can see from Tables VI and VII that the overall elec- 
tromagnetic field structure far away from the vertical conductor 
approaches that of TEM wave. 
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VERTICAL E-FIELD PRODUCED BY A 1-m LONG 
CURRENT DIPOLE 


In order to investigate the dependencies of scattered current 
magnitude on the thickness of vertical conductor and on the 
pulse width, we calculated vertical electric fields produced by a 
1-m long vertical current dipole shown in Fig. 19. This current 
dipole is energized by a single cycle, sinusoidal current pulse, 
which has a negative initial half cycle, a magnitude of | kA, 
and a total duration 7 = 200 ns (f = 5 MHz) or 100 ns (f = 
10 MHz). This current pulse is a rough approximation to the 
bipolar scattered current (Fig. 9). 


Upward E-fields 


Single-cycle 
” sinusoidal current 


Fig. 19. A 1m long vertical current dipole, on which a single cycle sinusoidal 
current pulse flows. 
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Fig. 20. Waveforms of vertical electric field produced by a single cycle sinu- 
soidal current pulse on a | m long vertical dipole: (a) f = 5 MHz (total duration 
T = 200 ns). (b) f = 10 MHz (7 = 100 ns). Downward electric field is defined 
as positive. 


Fig. 20(a) shows waveforms of upward vertical electric field 
at different horizontal distances from the dipole axis. It is clear 
from Fig. 20(a) that the magnitude of upward vertical electric 
field produced by the dipole decreases rapidly with increasing 
the horizontal distance. 

Fig. 20(b) is similar to Fig. 20(a) but for 7 = 100 ns (f = 
10MHz). It is clear from Fig. 20(b) that the magnitude 
of downward vertical electric field produced by the dipole 
decreases as the pulse becomes narrower. Within few meters 
from the 1-m current dipole, the static component of electric 
field, which is related to time integral of current, is dominant. 
Therefore, the electric field produced by the dipole at such 
horizontal distances decreases as the pulse gets narrower while 
its magnitude is the same. 
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